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Conversion  factors  for  U.S,  customary  to  metric  (SI)  units  of  measurement 


To  Convert  From 

To 

Multiply 

angstrom 

meters  (m) 

1.000  000XE-10 

atmosphere  (normal) 

kilo  pascal  (kPa) 

1.013  25  XE+2 

bar 

kilo  pascal  (kPa) 

1.000  000  XE+2 

barn 

meter2  (m2) 

1.000  000XE-28 

British  Thermal  unit  (thermochemical) 

joule  (J) 

1.054  350XE+3 

calorie  (thermochemical) 

joule  (J) 

4.184  000 

cal  (thermochemlcal)/cm2 

mega  Joule/m2(MJ  /m2) 

4.184  000XE-2 

curie 

giga  becquerel  (GBq)* 

3.700  000XE+1 

degree  (angle) 

radian  (rad) 

1.745  329  X  E-2 

degree  Fahrenheit 

degree  kelvln  (K) 

tK=(t°f  +  459.671/1.8 

electron  volt 

joule  (J) 

1.602  19  X  E- 19 

erg 

Joule  (J) 

1 .000  000  X  E-7 

erg/ second 

watt  (W) 

1.000  000  X  E-7 

foot 

meter  (m) 

3.048  000  X  E-l 

foot-pound-force 

joule  (J) 

1.355  818 

gallon  (U.S.  liquid) 

meter3  (m3) 

3.785  412  X  E-3 

inch 

meter  (m) 

2.540  000  X  E-2 

Jerk 

joule  (J) 

1.000  000XE+9 

Joule /kilogram  (J/Kg)  (radiation  dose 
absorbed) 

Gray  (Gy) 

1.000  000 

kllotons 

terajoules 

4.183 

kip  (1000  lbf) 

newton  (N) 

4  448  222  X  E^  3 

kip/inch2  (ksl) 

kilo  pascal  (kPa) 

6.894  757  X  E+3 

ktap 

newton-second/m2  (N-s/m2) 

1.000  000  XE+2 

micron 

meter  (m) 

1.000  000  X  E-6 

mil 

meter  (m) 

2.540  000  X  E-5 

mile  (international) 

meter  (m) 

1.609  344  X  E+3 

ounce 

kilogram  (kg) 

2.834  952  X  E-2 

pound-force  (lbf  avoirdupois) 

newton  (N) 

4.448  222 

pound-force  inch 

newton-meter  (N-m) 

1.129  848  X  E-l 

pound-force/lnch 

newton/meter  (N/m) 

1.751  268  X  E+2 

pound-force/ foot2 

kilo  pascal  (kPa) 

4.788  026  X  E-2 

pound-force/inch2  (psi) 

kilo  pascal  (kPa) 

6.894  757 

pound-mass  (ibm  avolidupois) 

kilogram  (kg) 

4.535  924  X  E-l 

pound-mass-foot2  (moment  of  inertia) 

kilogram-meter2  (kg  m2) 

4.214  C*i  1  X  E-2 

pound-mass/ foot3 

kilogram /meter3  (kg/m3) 

1.601  846  X  E+  1 

rad  (radiation  dose  absorbed) 

Gray  (Gy)** 

1.000  000  X  E-2 

roentgen 

coulomb /kilogram  (C/kg) 

2.579  760  X  E-4 

shake 

second  (s) 

1 .000  000  X  E-8 

slug 

kilogram  (kg) 

1.459  390  XE+1 

torr  (mm  Hg,  0°C) 

kilo  pascal  (kPa) 

1.333  22  X  E-l 
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••The  Gray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 
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SECTION  1 


INTRODUCTION 


The  amount  of  dust  and  debris  which 
can  be  scoured  and  lofted  into  a  nuclear 
cloud  depends  strongly  upon  surface 
conditions  in  the  target  area.  Regional 
and  seasonal  variations  in  land  use, 
vegetative  cover,  soil  moisture,  and  soil 
texture  all  influence  the  sweepup  mass 
source  and  lead  to  wide  differences  in 
the  mass  injection.  The  dependencies 
are  poorly  understood  and  are  not  at 
present  accounted  for  in  any  DNA 
sweepup  model.  As  a  result,  broad  un¬ 
certainties  are  introduced  any  time 
sweepup  models  are  applied  to  target 
areas  outside  the  region  for  which  they 
were  originally  developed  (namely,  dry. 
sandy  deserts).  This  is  a  serious  prob¬ 
lem,  especially  considering  the  geograph¬ 
ic  variety  of  strategic  targets  areas  found 
in  the  Soviet  Union  (e.g.,  subarctic  taiga 
forests,  black  soil  farmland,  wooded  river 
valleys,  and  semiarid  steppes).  Sweepup 
from  those  regions  can  hardly  be  ex¬ 
pected  to  resemble  sweepup  from  a  des¬ 
ert;  but  currently,  it  is  impossible  to 
distinguish  the  differences. 

In  this  report,  we  calculate  sweepup 
masses  (without  accounting  for  material 
fallback  to  the  surface)  for  a  wide  range 
of  surface  conditions.  Our  approach  has 
two  main  components:  (1)  A  surface 
classification  scheme  based  on  mea¬ 


sured  threshold  shear  velocities,  and  (2) 
a  theoretical  mass  flux  model  appropri¬ 
ate  for  the  high  surface  wind  speeds  as¬ 
sociated  with  nuclear  blast  waves.  These 
two  components  are  combined  such  that 
the  vertical  dust  flux  has  an  explicit  de¬ 
pendence  upon  the  threshold  shear  ve¬ 
locity.  Time-  and  range-dependent 
nuclear  blast  winds  are  prescribed  using 
analytic  positive  phase  dynamic  pres¬ 
sure  formulae.  Although  our  model  is 
highly  idealized  (it  cannot,  for  example, 
account  for  precursor  flows  or  uneven 
terrain),  the  results  nevertheless  demon¬ 
strate  the  seriousness  of  ignoring  region¬ 
al  and  seasonal  variations  in  target  area 
conditions.  Moreover,  they  indicate  a 
clear  direction  for  more  detailed  theoreti¬ 
cal  and  experimental  work. 

The  report  is  divided  into  six  sections.  In 
Sec.  2,  we  describe  the  surface  classifi¬ 
cation  scheme  and  discuss  how  thresh¬ 
old  velocities  can  vary  with  soil  type  and 
season.  Next,  in  Sec.  3,  we  describe  the 
sweepup  model.  This  model  is  then 
applied  in  Sec.  4  to  calculate  the  total 
sweepup  mass  and  the  "effective  scour¬ 
ing  radius”  for  various  surfaces,  heights 
of  burst,  and  yields.  The  implications  of 
these  results  are  discussed  in  Sec.  5, 
and  conclusions  presented  in  Sec.  6. 
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SECTION  2 


CLASSIFICATION  OF  SURFACE  TYPES 


The  blast  winds  generated  by  low  alti¬ 
tude  nuclear  detonations  are  strong 
enough  to  scour  material  from  virtually 
any  natural  surface.  Dirt,  pebbles,  bush¬ 
es — even  rocks  and  trees — can  be  carried 
away  when  subjected  to  winds  In  excess 
of  several  hundred  meters  per  second. 
Yet,  while  It  Is  clear  some  sweep  up  will 
probably  always  occur  (especially  at 
ground  ranges  where  dynamic  pressui  es 
maximize),  it  is  equally  clear  that  the 
amount  of  material  swept  up  is  not  the 
same  for  all  surfaces.  Vegetated  land,  for 
example,  yields  less  sweepup  mass  than 
barren  ground;  and  moist  (or  frozen)  soil 
less  than  dry.  The  key  Issue  Is  the  ability 
of  a  surface  to  resist  the  shear  stresses 
exerted  by  a  nuclear  blast  wave  and  Its 
accompanying  winds.  Any  factor  which 
either  (a)  lessens  the  stress  applied  di¬ 
rectly  to  the  soil  (such  as  a  plant  cover), 
or  (b)  Increases  the  cohesiveness  of  par¬ 
ticulates  lying  on  the  surface  (such  as  a 
high  moisture  content)  tends  to  lessen 
the  susceptibility  of  the  soil  to  sweepup 
and  thus  reduce  the  net  mass  lofted.  An 
explicit  quantification  of  all  those  factors 
has  not  been  developed;  however,  the 
threshold  shear  velocity  can  at  least 
roughly  categorize  their  net  effect  for  a 
given  surface. 

2.1  THRESHOLD  SHEAR  VELOCITY. 

The  threshold  shear  velocity  is  a  mea¬ 
sure  of  the  minimum  shear  stress  re¬ 
quired  to  initiate  particle  motion  along 
the  ground.  It  is  experimentally  deter¬ 
mined  in  the  field  using  a  portable,  op¬ 
en-bottomed  wind  tunnel  with  a 


regulated  airflow  |e.g.,  Gillette,  1978J.  As 
the  wind  speed  in  the  tunnel  is  in¬ 
creased  from  zero,  the  shear  velocity 


(1) 


is  monitored  until  saltation  (indicated  by 
particles  rolling  or  bouncing  along  the 
surface)  begins. 1  Here  r,  is  the  shear 

stress,  p  is  the  air  density,  k  *  0.4  is  the 
von  Karmann  coefficient,  and  U  is  the 
horizontal  wind  speed  at  some  height  Z 
above  the  ground.  The  value  of  U,  at 

which  saltation  first  occurs  defines  the 
threshold  shear  velocity  U,([) . 


The  threshold  shear  velocity  U  h  can 

also  be  used  to  define  the  point  at  which 
dust  enters  into  suspension.  This  regime, 
defined  by  the  condition  (Owen,  1964) 

-J71-  £  10  ,  (2) 

U*th 

begins  when  the  turbulent  drag  force 
upon  individual  particles  first  exceeds 
their  gravitational  fall  speed.  It  is  at  this 
point  that  the  particles  become  lofted 
and  can  be  carried  to  high  altitude.  Sus¬ 
pension  is  the  mechanism  which  creates 
sand  and  dust  storms  in  arid  regions.  It 
is  also  responsible  for  the  dense  sweep¬ 
up  layers  generated  by  nuclear  air- 
bursts.  Equation  (2)  therefore  forms  a 
base  for  the  sweepup  model  described  in 
Sec.  3. 


1.  The  relation  between  shear  stress  and  vertical  wind  shear  given  by  Eq.  (1)  is  strictly  valid  only  when 
the  atmospheric  surface  layer  is  statically  neutral  (well  mixed).  Corrections  to  this  formula  for  other 
stability  profiles  are  available  (see  Buslnger,  1973]. 
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2.2  OBSERVED  VARIATIONS. 

The  concept  of  using  U  h  as  an  index  for 

sweepup  potential  was  originally  devel¬ 
oped  for  application  to  "normal"  (non¬ 
nuclear)  conditions  such  as  dust  erosion 
from  drought-stricken  agricultural  land 
and  deserts.  Thus  the  only  systematic 
measurements  made  to  date  have  been 
over  dry,  barren,  or  sparsely  vegetated 
ground,  since  such  surfaces  are  most 
likely  to  be  subject  to  significant  low- 
wind  speed  sweepup.  Nevertheless,  even 
for  this  restricted  range  of  surface  types, 
wide  variations  in  U,((i  have  been  found. 

Table  1  summarizes  some  of  these  varia¬ 
tions.  In  general,  soils  which  are  loose 
(c.g.,  plowed  or  otherwise  disturbed)  or 
sandy  are  the  most  easily  eroded  (and 
thus  have  the  lowest  U  h );  loamy  soils, 

soils  with  a  large  number  jf  large  ele¬ 
ments  such  as  pebbles  and  those  with  a 
crusted  or  cloddy  surface  are  the  least 
erodibie.  In  most  cases,  mineral  content 
also  influences  soil  erosion  potential. 
High  concentrations  of  exchangeable  so¬ 
dium  or  calcium  carbonate  (CaC03),  for 


example,  tend  to  Increase  Ut(h  in  clay 

soils  by  strengthening  their  crustal 
strength  (Gillette,  1982). 

All  of  the  soils  listed  in  Table  1  were  dry. 
Technically,  this  means  that  the  water 
content  in  the  first  few  inches  below  the 
surface  was  at  or  below  the  wilting  point 
for  most  crops.  (The  wilting  point  corre¬ 
sponds  to  a  water  to  soil  volume  ratio  of 
0.05  for  sandy  loams,  0.08  for  loam,  and 
0.17  for  clay  [Alderfer,  1977].)  It  has 
been  found,  however,  that  once  the  wa¬ 
ter  content  of  a  soil  begins  to  exceed  the 
wilting  point  (due  to  rainfall  or  irriga¬ 
tion).  U,th  Increases  substantially  above 

the  values  listed  in  Table  1 ,  although  it 
is  not  clear  how  much  (Chepil,  1956]. 
Saturated  or  frozen  soils  are  likely  to  be 
even  more  stable,  but,  again,  no  mea¬ 
surements  are  available. 

Precipitation  can  modify  the  threshold 
velocity  in  other,  although  less  obvious, 
ways.  Sandy  soils  which  are  initially  dry 
and  cloddy,  for  example,  often  become 
more  erodable  several  days  to  weeks  af¬ 
ter  a  heavy  rainfall.  This  is  due  to  the 


Table  1 .  Mean  threshold  velocities  for  unvegetated  soils. 


Threshold 
Shear  Velocity 
(m/s) 

Percent  Mass 
in  Particles 
>  1  mm  diameter 

Modulus  of 
Crustal  Rupture 
(bars) 

Soil  Type 

Loose 

Cloddy 

Crusted 

Loose  Cloddy 

Sand 

0.28 

0.75 

0.66 

3 

60 

0.03 

Sandy  loam 

0.29 

1.05 

2.90 

30 

64 

0.42 

Loamy  sand 

0.34 

0.85 

1.03 

26 

47 

0.5 

Clay 

0.54 

>1.50 

>2.00 

42 

94 

0.75 

Silty  Clay 

0.56 

.... 

.... 

6 

-- 

.... 

Silty  clay 

0.64 

.... 

>1.50 

18 

-- 

.... 

loam 

Clay  loam 

0.68 

>1.09 

1.20 

28 

81 

0.38 

Loam 

0.78 

>1.50 

>1.50 

49 

89 

0.66 

Silt  loam 

1.08 

>2.00 

>1.50 

77 

89 

0.8 

Source:  Gillette  [1988]. 
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tendency  for  rainfall  to  “melt"  (disaggre¬ 
gate)  the  clods  lying  on  the  surface. 

Once  the  soil  dries,  it  returns  to  a 
smooth,  loose,  and  thus  more  easily 
eroded  state.  On  the  other  hand,  if  the 
soil  is  initially  loose,  wetting  and  drying 
may  ultimately  lead  to  the  formation  of  a 
crust,  thereby  increasing  U,(h.  Timing  is 

an  important  issue  here:  If  the  surface  is 
subjected  to  strong  winds  very  soon  after 
a  rainfall  (or  snowmelt),  significant  dry¬ 
ing  of  the  upper  few  millimeters  may  oc¬ 
cur  before  the  crust  can  form.  Sweepup 
can  then  follow,  even  though  the  soil  a 
centimeter  or  two  below  the  surface  is 
still  water-soaked.  This  effect  has  been 
observed  many  times  in  loam  and  clay 
loam  soils  [Gillette,  1988).  In  a  nuclear 
environment,  thermal  irradiation  of  the 
surface  may  also  dry  the  upper  soil  hori¬ 
zon,  possibly  negating  the  stabilizing  in¬ 
fluence  of  soil  moisture.  (This  influence 
could  however  be  mitigated  by  a  smoke, 
dust,  or  steam  layer  just  above  the  sur¬ 
face  which  would  tend  to  shield  the 
ground.) 

Very  few  measurements  have  been  made 
of  the  effects  of  live  vegetation  upon 
sweepup.  It  has  been  found,  howeve ', 
that  any  sort  of  grass  or  small  grain  crop 
cover  (e.g.,  wheat,  barley,  rye)  is  usually 
sufficient  to  preclude  sweepup  for  all  but 
the  strongest  natural  winds 

(ll,  >  2  m/s)  [Gillette,  1 988). 2  Even  a 

sparse  vegetation  cover  with  a  fractional 
lateral  cover  (frontal  silhouette  area  di¬ 
vided  by  ground  area)  greater  about  2  to 
4  percent  can  reduce  the  shear  stress 
upon  the  surface  to  near  zero  under  nat¬ 
ural  wind  conditions  [Marshall,  1971]. 
Nevertheless,  there  presumably  exists  an 


upper  bound  on  U,  above  which  sweepup 

commences.  This  critical  value  likely  de¬ 
pends  not  only  on  the  type  of  plant  and 
tfs  stand  density  (plants  per  unit  surface 
•rea),  but  also  on  the  plant  height,  its 
”  iter  content,  stem  thickness,  and  aero¬ 
dynamic  cross  section.  However,  for  very 
dense  plant  covers  such  as  small  grains 
or  grasses,  we  expect  that  the  plants 
would  either  be  lodged  (i.e.,  knocked 
over),  broken,  or  uprooted  long  before 
this  critical  shear  stress  Is  reached.  The 
surface  would  no  longer  be  as  well  pro¬ 
tected  and  the  critical  shear  velocity 
would  drop  to  a  value  more  representa¬ 
tive  of  bare  soils.  (Similar  modification 
would  result  if  the  plant  cover  were 
burned  due  to  pre-shock  thermal  irradi¬ 
ation.)  Unfortunately,  no  experimental 
data  exists  which  relates  high  speed 
sweepup  to  shear  stress  over  vegetated 
ground,  although  some  measurements  of 
the  critical  bending  moment  for  wheat 
and  barley  stem  breakage  have  been 
made  [Oda,  Suzuki,  and  Odagawa, 

1966],  Lacking  such  data,  we  arbitrarily 
assume  that  a  dense  plant  cover  results 
in  a  threshold  shear  velocity  of  at  least  3 
m/s.  Clearly,  experiments  must  be  per¬ 
formed  to  refine  this  value. 

2.3  SEASONAL  VARIATIONS. 

All  of  the  effects  discussed  above— the 
influence  of  vegetation  cover,  soil  mois¬ 
ture,  land  use,  and  soil  texture — imply 
that  the  potential  for  sweepup  (and 
thus  U,th )  must  vary  seasonally.  Figure 

1 ,  derived  from  a  database  used  to  pre¬ 
dict  soil  erosion  in  U.S.  agricultural 
areas  [Gillette  and  Passi,  1988],  illus¬ 
trates  the  expected  variation  U>th  for 


2.  One  exception  occurs  when  a  planted  field  lies  Immediately  downwind  from  a  barren,  sandy  (low 
U.lh)  surface  from  which  sweepup  is  occurring.  Material  from  the  upwind  source  can  then  spread 
into  the  field,  effectively  “sandblasting’*  the  plants  and  the  underlying  ground  surface,  thus  initiat¬ 
ing  sweepup  from  where  it  would  otherwise  not  be  expected  (Gillette,  1988].  Obviously,  such  down 
wind  sandblasting  effect  could  also  be  important  in  target  areas  containing  a  mixture  of  vegetated 
and  bare  [e.g..  fallow)  surfaces. 


4 


Figure  1.  Threshold  shear  velocity  as  function  of  season  for  U.S,  agricultural  area — 
loam  soil. 
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loam  soil  as  function  of  time  of  year  for 
two  land  use  types  (fallow  land  and  a 
spring-planted  crop)  and  for  two  “clima¬ 
tic”  classes  (corresponding  to  a  pro¬ 
longed  moist  period  and  a  prolonged 
drought,  respectively).  A  similar  plot  for 
sandy  soil  is  shown  in  Fig.  2.  Obviously, 
the  most  important  influence  is  the  pres¬ 
ence  and  maturity  of  the  crop  (in  this 
case,  assumed  to  be  a  small  grain  cere¬ 
al),  which  greatly  increases  the  stability 
of  the  soil  from  shortly  after  germination 
in  the  spring  until  harvest  in  late  sum¬ 
mer.  On  the  other  hand,  in  spring,  when 
the  ground  is  plowed,  crusts  and  clods 
are  broken,  the  threshold  velocity  is  re¬ 
duced,  and  the  soil  becomes  more  vul¬ 
nerable  to  wind  erosion.  After  harvest, 
threshold  velocities  remain  moderately 
high  due  to  the  common  practice  of  leav¬ 
ing  a  vegetative  residue  or  stubble 
through  the  winter.  In  drought  years, 
however,  Insufficient  moisture  is  avail¬ 
able  to  support  this  residue  and  clodding 
and  crusting  of  the  surface  is  inhibited. 
Threshold  velocities  therefore  remain 
low. 

2.4  NUCLEAR  INFLUENCES. 

Threshold  velocity  measurements  to  date 
have  focused  only  on  conditions  likely  to 
be  encountered  under  normal  condi¬ 
tions.  In  a  nuclear  environment,  howev¬ 
er,  several  other  factors  may  prove  to 


have  a  crucial  influence  upon  sweepup. 
For  example,  thermal  radiation  incident 
upon  the  surface  prior  to  arrival  of  the 
blast  wave  can  heat  the  soil  to  a  point 
where  bound  water  is  explosively  re¬ 
leased,  thus  causing  “popcorning”  of 
particles  from  the  surface  IVersteegen, 
Rault,  and  Hillendahl,  1989].  This  effect, 
which  is  most  prevalent  for  clay  soils, 
may  disrupt  surface  crusts  and  lower 
the  threshold  velocity  of  the  soil.  On  the 
other  hand,  thermal  irradiation  can  also 
melt  or  "glaze"  the  surface,  thereby  sta¬ 
bilizing  the  soil  by  leaving  a  thin  but 
presumably  strong  crust.  Sandy  sur¬ 
faces  are  probably  the  most  susceptible 
to  glazing,  although  it  has  only  rarely 
been  observed,  most  notably  following 
the  TRINITY  test  in  1945  [P.  Verstecgen. 
personal  communication,  1989], 

Blast  effects  can  also  modify  the  dust 
producing  potential  of  the  surface.  Air- 
blast  loading  of  subsurface  air  pores  can 
cause  spalling  and  disruption  of  an  un¬ 
broken  surface  with  the  passage  of  the 
blast  wave.  The  blast  can  also  uproot 
trees  and  bushes,  which  would  break 
the  surface  and  reduce  the  threshold  ve¬ 
locity,  as  well  as  increasing  the  shear 
stress  upon  the  ground.  Models  relating 
blast  effects  to  tree  blowdown  do  exist 
[Morris,  1973];  however,  no  quantifica¬ 
tion  of  its  effect  upon  soil  lofting  is  yet 
available. 
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SECTION  3 


SWEEPUP  MODEL 


In  the  previous  section,  we  reviewed 
many  of  the  links  between  surface  condi¬ 
tions  and  the  potential  for  dust  lofting.  A 
general  model  would  account  for  soil 
type,  texture,  moisture,  temperature  and 
vegetation.  In  this  first  analysis  we 
roughly  account  for  all  parameters  using 
a  single  quantitative  measure — the 
threshold  shear  velocity.  The  sweepup 
model  we  develop  is  clearly  an  approxi¬ 
mation,  but  nonetheless  accounts  for 
real  soils  in  calculating  the  amount  of 
dust  leaving  the  surface. 

3.1  BACKGROUND. 

Sweepup  models  relate  the  vertical  mass 
flux  F  of  dust  particles  leaving  the  sur¬ 
face  to  the  shear  stress  t*  =  p  U,  2  ex¬ 
erted  on  the  ground  by  a  sheared 
(rotational)  wind  field.  The  general  form 
is 


f*«  (u.-uJ,i  131 

where  U,  >  U  h  and  a  and  ft  are  empir¬ 
ical  coefficients.  At  low  wind  speeds,  di¬ 
mensional  analysis  suggests  that  F 
should  vary  with  the  kinetic  energy  de¬ 
livered  to  the  surface  [Gillette,  1980). 
Thus  we  expect  ft  ~  3;  this  dependence 

has  indeed  been  borne  out  by  observa¬ 
tions  of  natural  sweepup  from  desert 
surfaces  [Shinn  et  al.,  1976;  Westphal, 
Toon,  and  Carlson,  1987).  However,  as 
wind  speeds  rise  and  dust  concentra¬ 
tions  increase,  the  lofted  particles  begin 
to  act  as  a  momentum  sink  on  the  air¬ 
flow.  Thus  the  shear  on  the  surface  is 
reduced,  the  lofting  ability  of  the  wind 
diminished,  and  F  is  no  longer  propor¬ 
tional  to  U,3.  Unfortunately,  there  are 


very  few  measurements  of  F  in  this  re¬ 
gime,  especially  at  the  very  high  wind 
speeds  characteristic  of  nuclear  blast 
waves.  In  one  experiment  [Hartenbaum, 
1971],  a  value  ft  ~  1.0  was  found  over 
an  uncohesive  sand  (mean  particle  ra¬ 
dius  125//m;U,th  ~  0.35  m/s)  surface 

for  free  stream  wind  speeds  U  between 
34  and  1 15  m/s  (corresponding  to  shear 
velocities  In  the  range  5.0  <  U,  ^  18.6 

m/s).  This  result  cannot  be  easily  ex¬ 
plained  by  simple  dimensional  analyses. 
A  model  which  accounts  for  the  effect  of 
dust  loading  on  the  airstream  at  very 
high  wind  speeds  is  required. 

3.2  HIGH  SPEED  MASS  FLUX  MODEL. 

For  wind  speeds  in  the  regime 

U,  >  10  U  h ,  the  vertical  mass  flux  is 

related  to  the  surface  shear  stress  by 
(Mirels,  1984); 


F  =  P0U. 


2  l  n(l  +  B) 
Uo 


(4) 


where  po  and  Uo  are  the  density  and 
horizontal  wind  velocity  in  the  free 
stream  (i.e.,  at  the  top  of  the  dusty  sur¬ 
face  layer)  and  B  is  a  dimensionless 
blowing  parameter  which  expresses  the 
effect  of  transverse  particle  injection  on 
the  local  shear  stress  in  plane  parallel 
flow.  Assuming  that  the  mass  loading  ef¬ 
fectively  reduces  the  shear  stress  to  the 
threshold  value  required  to  maintain 
dust  lofting,  B  is  given  by  the  implicit  re¬ 
lation 


Jln(l  +  B)  =  B 


(5) 
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which,  in  the  suspension  regime 
(U,  >  10  Unh),  can  be  approximated 


An(l  +  B)  -  3.92 


0.24 


(6) 


The  functional  dependence  of  B  on  the 
shear  stress  ratio  r,/r  h  given  by  Eqs.  (5) 

and  (6)  is  illustrated  in  Fig.  3. 


To  prescribe  Uo  as  a  function  of  time, 
ground  range,  yield,  and  height  of  burst, 
we  use  the  ideal  airblast  approximations 
given  by  Brode  (1987).  These  approxima¬ 
tions  are  valid  only  during  the  positive 
overpressure  phase  of  the  blast;  we 
therefore  assume  all  the  sweepup  occurs 
during  the  positive  phase.  The  free- 
stream  velocities  are  then  used  to  com¬ 
pute  the  shear  velocity  U,  according  to 


the  “rough  plate”  formula  (Schlichting, 
19581: 

[2.87  +  1.58  log(x/Ks)]-1,25  , 

(7) 

where  x  is  the  distance  behind  the  shock 
wave,  Ks  -  0.05  m  is  the  roughness 
height,  and  x  >  100  Ks.  Finally,  substi¬ 
tuting  Eqs.  (6)  and  (7)  into  Eq.  (4),  we  ob¬ 
tain  (see  Eq.  (8)  below! 

The  sweepup  mass  flux  is  therefore  in¬ 
versely  (but  weakly)  dependent  upon 
threshold  velocity.  This  expression  com¬ 
pares  well  with  the  experimental  result  F 
a  Uo1  144  derived  by  Hartenbaum  (1971) 
and  subsequently  used  as  the  basis  for 
many  nuclear  sweepup  models 
ISchlamp,  Schuckman,  and  Rosenblatt, 
1982;  Bacon,  Dunn,  and  Sarma,  1988J. 
However,  unlike  previous  models,  the  de¬ 
pendence  of  sweepup  on  threshold  veloc¬ 
ity  (and  hence  on  surface  type)  is  explicit 
in  Eq.  (8). 


1.80 


F  « 


Po 


U°;J.4  [2.87+  1.58  log(x/Ks)] 


2.80 


U0 


1.24 


th 


(8) 


9 


Blowing  parameter,  B 


Figure  3.  Blowing  parameter  as  function  of  normalized  shear  stress. 
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SECTION  4 


SWEEPUP  CALCULATIONS 


We  now  use  the  mass  flux  model  and 
ideal  blast  wave  driver  to  compute  the 
variation  in  net  sweepup  as  a  function  of 
threshold  shear  velocity.  Two  quantities 
are  of  interest:  ( 1 )  The  integrated  sweep¬ 
up  mass  M;  and  (2)  the  effective  sweepup 
radius  Re.  defined  as  the  the  maximum 
range  (from  ground  zero)  at  which 
sweepup  can  occur  for  a  given  yield. 

HOB.  and  surface  type. 

4.1  METHODOLOGY. 

We  treat  the  surface  as  flat,  homoge¬ 
neous  (i.e.,  no  spatial  variations  in 
threshold  velocity),  and  thermally  ideal. 
Assuming  cylindrical  symmetry  about 
the  burst  point,  the  incremental  dust 
mass  AM(  lofted  from  a  ring  of  width  AR 
and  radius  R,  is 

fT  +  Du 

AM^lTRiARj^  Fjdt  (9) 

where  T  is  the  blast  wave  time  of  arrival 
at  Rt,  Du  is  the  dynamic  pressure  posi¬ 
tive  phase  duration,  and  Ftis  the  local 
sweepup  mass  flux  (Eq.  8).  We  set  ARi . 
60  m  and  use  empirical  formulae  to  de¬ 
termine  T  and  Du  as  functions  of  height 
of  burst,  yield,  and  range  [Brode,  1987]. 
To  simulate  the  low  windspeed  cutoff  for 
dust  suspension,  we  further  assume  Ft  = 
0  for  U.  <  10U,th. 

Equation  (9)  is  integrated  numerically 
using  the  trapezoidal  rule  and  a  timestep 
At  =  Du/ 100  for  i  =  1  to  i  =  IMAX,  where 
IMAX  corresponds  to  the  first  ring  where 


no  sweepup  occurs  (AMjmax  =  0)-  Thus 
Rimax-i  =  Re  defines  the  effective  sweep¬ 
up  radius.  The  total  sweepup  mass  M  is 
then  found  by  simply  summing  over  the 
rings: 

IMAX 

M=  X  AMj  .  (10) 

i=i 

Note  that  M  is  a  measure  of  the  initial 
sweepup:  it  does  not  account  for  mass 
fallback  to  the  surface,  nor  does  it  in¬ 
clude  mass  lofted  during  the  negative 
phase. 

4.2  VARIATIONS  IN  SWEEPUP  MASS. 

We  have  computed  the  variation  in  total 
sweepup  mass  M  for  threshold  shear  ve¬ 
locities  from  0.20  to  4.0  m/s.  scaled 
heights  of  burst  (SHOBs)  from  20  to 
1000  ft/KTl/3,  and  yields  (W)  from  500 
to  1000  KT.3  Figure  4  shows  the  depen¬ 
dence  of  M  on  U,(h  for  a  500  KT  detona¬ 
tion  at  three  burst  altitudes:  397  ft 
(SHOB  =  50  ft/KTl/3),  1984  ft  (SHOB  = 
250  ft/KT,/3).  and  3969  ft  (SHOB  =  500 
ft/KT1  /3).  We  find  that  M  decays  expo¬ 
nentially  with  increasing  U,(h ,  and  is 

most  sensitive  to  variations  in  threshold 
velocity  at  low  U,(h  — i.e.,  over  dry,  loose, 

unvegetated  soils.  In  this  regime,  small 
changes  in  surface  type  can  produce 
large  changes  in  sweepup.  For  example, 
we  find  that  a  50  ft  SHOB  detonation 
over  loose  sand  ( U,th  =  0.28  m/s)  raises 

over  two  times  as  much  mass  as  an 


3.  The  burst  altitude  SHOB  =  20  ft/KT13  roughly  corresponds  to  the  minimum  for  non-cratering  air- 
bursts  (Rosenblatt,  1981], 
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50  ft  SHOB 
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Threshold  shear  velocity,  U.th  (m/s) 


identical  burst  over  loose  loam  ( U  h  = 

0.78  m/s).  If  the  burst  occurs  at  an  even 
higher  altitude  (e.g.,  500  ft  SHOB).  the 
difference  is  more  pronounced.  For  high¬ 
ly  non-erodible  soils  ( Ut(h  >  2.0)  or  soils 

which  are  densely  vegetated 
(U,(h  S  3.0),  sweepup  is  strongly  sup¬ 
pressed,  especially  at  higher  burst  alti¬ 
tudes.  Indeed,  we  find  that  when  a  500  ft 
SHOB  detonation  occurs  over  a  surface 
with  U  h  >  3.5  m/s  (characteristic  of  an 

early  to  mid-summer  wheat  crop — see 
Figs.  1  and  2),  no  sweepup  occurs  at  all. 

The  strong  dependence  of  sweepup  mass 
on  burst  altitude  is  illustrated  in  Fig.  5. 
Here  we  have  plotted  M  as  a  function  of 
SHOB  for  a  500  KT  burst  over  four  sur¬ 
face  types  ranging  from  loose  sand  to  a 
desert  pebble  pavement  (mass  equivalent 
modal  pebble  diameter  1 .0  cm)  similar  to 
that  found  in  Yucca  Valley  at  the  Nevada 
Test  Site  (NTS).  We  find  that  the  depen¬ 
dence  of  M  on  SHOB  is  highly  nonlinear, 
with  the  largest  variations  occurring  over 
the  most  erodible  soils.  In  particular,  we 
find  that  sweepup  mass  at  first  decays 
rapidly  with  increasing  SHOB,  reaching 
a  relative  minimum  around  80  to  100  ft 
SHOB.  With  further  increases  in  burst 
altitude,  however,  the  sweepup  mass  be¬ 
gins  to  rise — indicative  of  the  build-up  in 
blast  dynamic  pressure  with  increasing 
SHOB  for  shocks  waves  in  the  Mach  re¬ 
flection  regime.  The  positive  dependence 
of  M  on  SHOB  continues  until  an  alti¬ 
tude  is  reached  at  which  the  Mach  wave 
begins  to  weaken— typically  between  300 
and  400  ft  SHOB.  For  bursts  above  this 
altitude,  sweepup  again  diminishes, 
eventually  dropping  to  zero.  Over  desert 
pavement,  the  cutoff  occurs  at  SHOB  = 
760  ft/KT1/3  (HOB  =  6032  ft);  over  loam, 
it  is  at  SHOB  =  935  ft/KT1/3(HOB  = 

7421  ft).  These  results  are  modified 
somewhat  over  non-ideal  surfaces  due  to 
enhanced  velocities  in  the  precursed 
wave.  In  particular,  M  is  increased  rela¬ 


tive  to  the  results  in  Fig.  5  for  all 
SHOB  <  650  ft/KT1/3, 

Sweepup  also  increases  with  yield.  Fig¬ 
ure  6  shows  the  variation  in  M  for  four 
surface  types  and  two  burst  altitudes  (50 
and  250  ft  SHOB)  for  yields  between  500 
KT  and  1  MT.  The  calculations  show  a 
linear  dependence  on  yield  in  this  range, 
with  sandy  surfaces  being  the  most  sen¬ 
sitive  to  changes  in  W  and  desert  pave¬ 
ment  the  least.  Thus  the  sweepup 
efficiency  (defined  as  soil  mass  lofted  per 
megaton  of  weapon  yield)  is  yield  depen¬ 
dent.  For  example,  a  500  KT  burst  at  50 
ft  SHOB  produces  sweepup  efficiencies 
ranging  from  0.12  Tg/MT  for  desert 
pavement  to  0.33  Tg/MT  for  sand.  Rais¬ 
ing  the  yield  to  1  MT,  however,  nearly 
doubles  these  values  to  between  0.23 
and  0.64  Tg/MT.  Slightly  lower  efficien¬ 
cies  result  at  all  yields  when  the  burst 
occurs  at  a  higher  altitude  (250  ft  SHOB) 
over  the  two  least  erodible  surfaces 
(loam  and  desert  pavement),  while  higher 
efficiencies  are  obtained  over  the  most 
erodible  ones.  This  is  consistent  with 
Fig.  5. 

4.3  VARIATIONS  IN  SWEEPUP 
RADIUS. 

The  variations  in  sweepup  mass  dis¬ 
cussed  above  can  be  attributed  to  a  com¬ 
bination  of  two  factors:  (1)  Variations  in 
positive-phase  dynamic  impulse  inte¬ 
grated  over  the  sweepup  area,  and  (2) 
variations  in  the  size  of  the  sweepup 
area  itself.  Sweepup  area  can  be  charac¬ 
terized  by  the  effective  sweepup  radius 
Re.  Figure  7  shows  the  dependence  of  RE 
on  threshold  shear  velocity  for  three 
burst  altitudes.  Comparing  these  results 
with  Fig.  4,  it  is  clear  that  the  sweepup 
mass  M  depends  strongly  upon  RE— both 
decay  exponentially  with  increasing  U  h, 

and  both  display  the  greatest  sensitivity 
to  threshold  velocity  over  bare,  loose  soil. 
(Not  all  the  dust  to  radius  RE  is  lofted  to 
stabilization;  but  all  is  raised  from  the 
surface  and  is  directed  toward  the 
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Figure  5.  Net  sweepup  mass  as  function  of  burst  altitude. 
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Threshold  shear  velocity,  UMh  (m/s) 


pedestal.)  Unlike  the  sweepup  mass, 
however,  there  is  no  “crossover”  with  in¬ 
creasing  threshold  velocity;  the  sweepup 
radius  for  a  500  ft  SHOB  burst  is  always 
larger  than  for  detonations  at  lower  alti¬ 
tudes  (at  least  until  the  cutoff  at  U*th  = 

3.5  m/s  is  reached).  This  is  as  we  would 
expect  for  a  Mach  reflected  shock.  It  sim¬ 
ply  indicates  that  as  bursts  occur  at 
greater  heights  above  the  surface,  more 
of  the  kinetic  energy  is  converted  into 
horizontal  dynamic  pressure.  The  rela¬ 
tionship  between  RE  on  SHOB  is  illus¬ 
trated  in  more  detail  in  Fig.  8.  Unlike  the 
sweepup  mass  M  (Fig.  5),  we  find  that  RE 
varies  linearly  with  SHOB. 

Finally,  in  Fig.  9.  we  show  the  variation 
in  Re  as  a  function  of  yield.  Again,  the 
dependence  is  linear,  with  the  greatest 
sensitivity  found  for  the  most  erodible 
soils.  However,  when  expressed  in  terms 

of  sweepup  area  [n  RE  2) ,  this  sensitivity 

disappears.  For  example,  from  Fig.  9,  we 
find  that  a  1-MT  burst  always  sweeps  an 
area  about  1 .6  times  larger  than  a 
500-KT  detonation  at  the  same  altitude, 
regardless  of  soil  type. 

4.4  NEVADA  TEST  SITE  SOIL. 

All  of  the  results  presented  so  far  have 
implicitly  assumed  that  the  soil  proper¬ 
ties  are  independent  of  depth.  Thus  we 
have  assumed  that  the  threshold  shear 
velocity  remains  constant  throughout 
the  period  sweepup  is  occurring.  For 
many  agricultural  soils,  which  have  a 
deep,  well-mixed  upper  layer  of  topsoil, 
this  is  an  appropriate  assumption.  But 
many  natural  soils  are  not  so  homoge¬ 
neously  structured.  The  Yucca  Valley  re¬ 
gion  of  the  NTS  is  an  example.  Soils  in 
most  of  this  area  are  representative  of 
“wind-stabilized"  alluvium  from  which 
most  of  the  erodible  (i.e.,  small)  surface 
elements  have  long  since  been  removed 
by  natural  wind  processes.  As  a  result. 


the  soil  has  a  two-layer  structure:  A  thin 
desert  pavement  veneer  consisting  of 
large  pebbles  (tens  of  millimeters  to  a 
centimeter  or  two  in  diameter),  overlying 
a  deeper  layer  of  finer  material — typically 
sand  or  sandy  loam  imbedded  with  grav¬ 
el.  This  desert  pavement  is  usually  suffi¬ 
cient  to  prevent  wind  erosion  of  the  fine 
particles  for  all  but  the  strongest  natural 
winds;  in  a  nuclear  environment,  howev¬ 
er,  it  is  easily  removed,  as  has  been  ob¬ 
served  following  bursts  over  Yucca  Valley 
[Lamar,  1962],  It  therefore  seems  likely 
that  nuclear  sweepup  over  desert  pave¬ 
ments  may  proceed  in  two  stages:  A  first 
stage  during  which  only  the  large  par¬ 
ticles  comprising  the  top  layer  layer  are 
removed,  followed  by  second,  more  vigor¬ 
ous  stage  during  which  the  now-exposed 
underlying  soil  is  scoured.  Thus,  more 
mass  may  ultimately  be  lofted  than  if 
only  the  heavy  top  layer  were  present. 

To  test  whether  this  two-stage  process 
can  significantly  increase  the  net  sweep¬ 
up  mass,  we  have  modified  our  model  to 
account  for  the  presence  of  desert  pave¬ 
ment  over  a  fine  soil.  We  assume  that 
the  surface  is  initially  covered  with  a 
single  layer  of  pebbles  with  mass  modal 
diameters  Dm  =  1  cm.  To  this  layer  we 
assign  a  threshold  shear  velocity  U,(h  = 

1.36  m/s.  This  value  is  consistent  with 
the  empirical  relation  for  dry,  disturbed 
desert  soils  (Gillette  et  al..  1980) 

UMh  =  0.43  +  0.93  Dm  (11) 

where  Dm  is  in  centimeters.  Assuming 
the  particles  are  spherical  and  have  a 
mass  density  2.0  x  103  kg/m3,  we  calcu¬ 
late  the  areal  density  of  the  pebble  layer 
to  be  10.5  Ap  kg/m2,  where  0  <  Ap  s  1.0 
is  the  area  fraction  of  ground  actually 
covered  by  pebbles.  (We  estimate  Ap  - 
0.3  for  Yucca  Valley.)  Setting U,(h  =  1.36 

m/s  in  Eq.  8  and  using  Eq.  9,  the  time- 
dependent  pebble  mass  removed  from 
each  concentric  ring  around  ground  zero 
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is  calculated.  If,  before  the  sweepup 
phase  ends,  we  find  the  total  mass  re¬ 
moved  has  exceeded  the  critical  density 
10.5  Ap  kg/m2,  we  assume  that  the 
pebble  layer  has  been  completely  re¬ 
moved  from  that  ring.  We  then  “instanta¬ 
neously"  reset  U  h  to  a  lower  value  to 

simulate  the  exposure  of  the  underlying 
fine  particles,  and  proceed  with  the  cal¬ 
culation  until  sweepup  ends. 

We  have  computed  the  net  sweepup 
mass  for  a  50  ft  SHOB  burst  at  two 
yields  (500  and  1000  KT)  over  three  dif¬ 
ferent  desert  pavement  configurations: 

( 1 )  A  single  pebble  veneer  over  highly 
erodible  sand  (U,(h  =  0.28  m/s),  (2)  a 

single  pebble  veneer  over  an  intermedi¬ 
ately  erodible  soil  ( U,lh  =  0.45  m/s),  and 


(3)  an  “infinitely  deep”  pebble  layer  for 
which  U  h  remains  fixed  at  1.36  m/s 

throughout  the  calculation.  The  results 
are  presented  in  Table  2.  We  find  desert 
pavement  to  be  surprisingly  more  stabi¬ 
lizing  than  expected.  Only  when  the 
pavement  overlies  loose  sand  and  is 
sparsely  distributed  (Ap  =  0.3)  do  we  find 
an  enhancement  in  sweepup  mass  in  ex¬ 
cess  of  10  percent.  This  is  in  spite  of  the 
fact  that  we  compute  complete  removal 
of  the  pebble  veneer  (for  An  =  0.3)  out  to 
approximately  425  ft/KT1'3  range  (~1.4 
km  for  a  1-MT  burst).  This  is  equivalent 
to  nearly  one-quarter  the  entire  sweepup 
area.  However,  the  most  significant 
scouring  of  the  underlying  soil  is  re¬ 
stricted  to  ranges  much  closer  to  ground 
zero;  thus  the  net  sweepup  mass  is  pro¬ 
portionately  reduced. 


Table  2.  Sensitivity  of  sweepup  mass  to  desert  pavement  cover  for  50  ft  SHOB  burst. 


Threshold 
Shear  Velocity 
of  Underlying 
Soil  (m/s) 

Yield.  W 
(KT) 

Areal 

Pebble 

Cover, 

Ap 

Sweepup  Mass,  M 
~  (kt) 

Pebbles 

Pebbles  Overlying 

Only  Fine  Soil 

Difference 

(%) 

0.28 

500 

0.3 

60.09 

65.77 

+9.42 

0.5 

60.09 

63.67 

+5.92 

1.0 

60.09 

62.31 

+3.66 

1000 

0.3 

117.20 

130.30 

+  11.18 

0.5 

117.20 

125.70 

+7.25 

1.0 

117.20 

122.40 

+4.44 

0.45 

500 

0.3 

60.09 

63.90 

+6.  34 

0.5 

60.09 

62.36 

+3.78 

1.0 

60.09 

60.10 

+0.02 

1000 

0.3 

117.20 

126.20 

+7.68 

0.5 

117.20 

122.70 

+4.69 

1.0 

117.20 

120.40 

+2.73 
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These  results  would  be  drastically  differ¬ 
ent  for  non-ideal  surfaces  when  thermal 
layers  and  precursors  are  accounted  for. 
The  resultant  elevated  dynamic  pressur¬ 
es  would  likely  remove  the  pebble  veneer 
much  sooner  and  to  a  greater  distance 
from  ground  zero.  The  net  sweepup 
mass  would  therefore  be  higher.  This  is 
partially  supported  by  the  observation  of 
Lamar  (1962),  who  found  near  complete 
removal  of  the  desert  pavement  out  to  at 
least  600  ft/KT1/3  from  ground  zero  fol¬ 
lowing  shot  SHASTA  (W  =  16.5  KT, 

SHOB  =  196  ft/KT1  /3)  in  Yucca  Valley. 
Separate  evidence,  however,  indicates 
that  this  shot  did  not  produce  a  strong 
precursor  [Liner  et  al..  19751. 

4.5  NORMALIZED  SWEEPUP  MASS. 

Figure  10  shows  the  variation  in  sweep¬ 
up  mass  as  a  function  of  threshold  ve¬ 
locity  normalized  by  that  calculated  for 
Yucca  Valley.  NTS.  To  compute  the  nor¬ 
malization  factor  Mnts,  we  adopted  the 
“two-layer"  model  described  above  with 
Ap  =  0.3  and  U,(h  =  0.45  m/s  for  the  un¬ 


derlying  soil;  the  results  are  indicated  in 
the  figure.  (Our  choice  U  h  =  0.45  m/s  is 

somewhat  arbitrary;  it  probably  underes¬ 
timates  the  dust  pickup  at  the  NTS  since 
some  of  the  desert  pavement  was  broken 
by  pre-shot  construction  and  traffic.) 
These  results  emphasize  that,  compared 
to  most  dry  barren  or  semi-barren  dry 
surfaces,  Yucca  Valley  soil  is  not  as 
readily  eroded.  This  is  due  to  the  stabi¬ 
lizing  effect  of  the  desert  pavement 
pebble  cover.  Most  other  bare  soils — for 
example,  agricultural  soils  during  a  dry 
winter  or  immediately  after  plowing  or 
harvest — actually  produces  more  sweep¬ 
up  than  the  Yucca  Valley  surface.  For 
example,  for  detonations  at  50  to  500  ft 
SHOB  over  bare  loam  surfaces  (  U,(h  = 

0.78  m/s),  between  1.3  and  2.2  times 
more  dust  would  be  raised.  On  the  other 
hand,  if  the  surface  were  heavily  vege¬ 
tated  (  U  h  >  3.5  m/s),  the  sweepup  is 

reduced  to  anywhere  from  zero  to  about 
50  percent  of  its  NTS  mass. 
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SECTION  5 


IMPLICATIONS  FOR  TARGET  AREAS 


Surfaces  in  target  areas  vary  widely  re¬ 
flecting  differences  in  soil  type,  topogra¬ 
phy,  moisture  levels,  and  vegetation 
cover.  Most  real  target  surfaces  do  not 
resemble  the  Nevada  desert,  nor  do  they 
remain  invariant  throughout  the  year. 
Small  differences  in  surface  properties 
can  lead  to  large  changes  in  dust  sweep- 
up. 

The  real  soil  sweepup  model  we  devel¬ 
oped  used  a  single  parameter,  the 
threshold  velocity,  to  account  for  soil 
type,  moisture,  and  vegetation  cover. 
Clearly  this  is  a  simplification,  but  one 
that  nonetheless  accounts  for  exper¬ 
imental  data  on  low-speed  scouring  of 
real  soils.  In  extending  this  analysis  to 
nuclear  sweepup,  we  recognize  that  sev¬ 
eral  potentially  important  physical  pro¬ 
cesses  are  treated  inadequately  or  not  at 
all.  Inclusion  of  (a)  non-ideal  shock  ef¬ 
fects,  (b)  blast  and  thermal  modification 
of  surface  properties,  (c)  uneven  terrain 
effects,  and  (d)  post-positive  phase  blast 
winds  would  no  doubt  alter  the  results 
somewhat.  We  anticipate  that  each  of 
these  processes  modifies  the  real  soil 
corrections  we  calculate,  nevertheless, 
the  conclusions  indicate: 

•  Targets  located  on  cropland  have  the 
greatest  dust  producing  potential, 
since  the  surface  tends  to  be  bare  at 
least  part  of  the  year— especially  dur¬ 
ing  the  spring  plowing  and  immedi¬ 
ately  after  harvest  in  the  autumn 
(Figs.  1  and  2).  Moreover,  some  por¬ 
tion  of  the  land  remains  bare  (i.e., 
fallow)  all  year  around.  For  example, 
in  U.S.  target  regions,  typically  one- 
fourth  the  area  devoted  to  crops  is 
fallow  at  any  given  time  [North  Dako¬ 
ta  Agricultural  Statistics  Service, 
1988], 


•  A  crop  or  grass  cover  may  be  suffi¬ 
cient  to  completely  suppress  sweepup 
from  a  500  KT  airburst  if  the  SHOB 
is  greater  than  500  ft/KT1/3  (Fig.  4). 
Thus  late  summer  is  the  least  vulner¬ 
able  time  of  year  for  sweepup:  this  is 
when  crop  and  grass  densities  are 
greatest. 

•  There  exists  an  optimum  burst  alti¬ 
tude  between  SHOB  =  80  to  100  ft/ 
KT1/3  at  which  the  sweepup  mass  is 
minimized  over  bare  soils.  Bursts 
detonated  below  this  altitude  or  be¬ 
tween  300  and  400  ft/KTl/3  tend  to 
maximize  dust  production  (Fig.  5). 

•  For  a  fixed  SHOB,  sweepup  mass 
tends  to  increase  with  increasing 
yield.  The  sensitivity  to  yield  is  great¬ 
est  for  the  most  erodible  soils  (Fig.  6). 

•  The  effective  sweepup  radius  in¬ 
creases  with  increasing  burst  altitude 
and  yield,  and  is  most  sensitive  to 
surface  conditions  over  highly  erod¬ 
ible  soil  (Figs.  7-9). 

•  Bare,  dry  agricultural  soils  (com¬ 
posed  primarily  of  loam)  produce  up 
to  220  percent  more  sweepup  mass 
than  a  desert  pavement  surface  typi¬ 
cal  of  the  NTS  (Fig.  10). 

Our  results  also  imply  that  sweepup 
may  not  be  uniform  even  within  a  single 
target  area.  Intercontinental  Ballistic 
Missile  silo  fields  are  quite  large,  typical¬ 
ly  covering  an  area  on  the  order  and  1 04 
km2  in  the  U.S.  and  somewhat  less  in 
the  Soviet  Union.  Surface  conditions 
within  an  area  so  large  are  not  necessar¬ 
ily  uniform.  A  burst  over  one  portion  of  a 
silo  field  could  yield  a  much  different 
sweepup  amount  and  scour  an  area  con¬ 
siderably  larger  or  smaller  than  a  burst 
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in  another  portion  of  the  field.  Fratricide 
probabilities  could  thus  be  highly  vari¬ 
able. 

Figure  1 1  depicts  a  U.S.  target  area. 
Aside  from  isolated  urban  areas  and  bo¬ 
dies  of  water,  most  of  the  land  shown 
here  is  either  planted  (primarily  with 
wheat,  oats,  and  barley)  or  is  tail-grass 
wildland.  Yet  the  soil  upon  which  this 
vegetation  grows  is  not  the  same  across 
the  entire  region;  roughly  half  the  silos 
are  located  on  loamy  soil  which  is  well- 
drained  and  hence  dry  (except  for  peri¬ 
ods  immediately  after  a  rainfall  or  in 
irrigated  zones).  Thus,  assuming  U  h  = 


0.78  m/s,  our  results  indicate  a  maxi¬ 
mum  sweepup  efficiency  of  about  0. 16 
Tg/MT  for  a  500-KT  low  altitude  burst, 
ranging  up  to  ~  0.30  Tg/MT  for  a  1-MT 
detonations.  In  addition,  we  calculate 
maximum  sweepup  radii  between  3  and 
4  km  from  the  burst  point — roughly  one- 
third  to  one-half  the  mean  spacing  be¬ 
tween  silos.  Smaller  sweepup  radii  and 
lower  efficiencies  are  to  be  expected  from 
strikes  against  those  silos  located  on 
moist  or  poorly  drained  soils,  or  follow¬ 
ing  attacks  during  less  vulnerable  times 
of  year  such  as  summer  or  when  the 
surface  is  snow-covered  or  frozen. 


Source:  Aandahl  [1982];  Omodt  et  al.  [1968], 


Figure  11.  Soils  in  the  Minot  AFB,  North  Dakota  ICBM  silo  area 
(dots  indicate  silo  locations). 
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SECTION  6 


CONCLUSIONS 


We  find  that  the  threshold  shear  velocity 
provides  a  useful  quantitative  measure 
of  the  susceptibility  of  a  surface  to  nu¬ 
clear  sweepup.  Wide  variations  in  this 
parameter  occur  naturally  with  seasonal 
and  regional  differences  in  vegetation 
cover,  soil  moisture,  soil  type,  and  land 
use.  We  find  that  such  variations  can 
lead  to  large  changes  in  the  initial 
sweepup  mass  raised  by  a  low  altitude 
airburst — particularly  over  surfaces 
which  are  dry  and  minimally  vegetated. 
When  vegetation  is  present,  the  surface 
is  strongly  stabilized;  in  fact,  a  dense 
agricultural  cover  can  in  some  cases 


completely  suppress  sweepup.  By  ex¬ 
trapolating  our  results  to  real  soil  condi¬ 
tions  in  target  regions,  we  expect  even 
larger  variations  in  sweepup  mass  could 
result  under  attack  conditions.  We  find 
that  sweepup  is  minimal  as  long  as  the 
vegetation  cover  remains  intact;  at  high 
overpressures  or  high  thermal  fluxes, 
the  stabilizing  influence  is  diminished. 
The  threshold  velocities  at  which  plant 
lodging,  breakage,  and  uprooting  occur 
are  not  yet  established,  but  it  is  never¬ 
theless  clear  that  sweepup  is  greatly  re¬ 
duced. 
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